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Priority Data 

This application is a divisional application of U.S. patent application Serial No. 
10/037,620, titled "Estimating A Plurality of Tones in an Input Signal" filed January 3, 
2002, whose inventor is Yong Rao, which is a continuation-in-part of U.S. patent 
10 application Serial No. 09/880,167 titled "System and Method for Estimating One or More 
Tones in an Input Signal" filed June 12, 2001, whose mventors are Yong Rao and Michael 
Cema. 

Field of the Invention 

15 The invention relates generally to the field of signal analysis, and more 

particularly, to detecting the fi-equency, amplitude and/or phase of a plurality of tones 
comprised within an input signal. 

Description of the Related Art 

20 In many applications, it is necessary or desirable to precisely locate multiple tones 

comprised in a signal. This need arises in many fields including telecommimications, 
radar, sonar, networking, measurement, and various other applications. Prior art 
techniques for detecting tones in a signal generally may not produce accurate results 
and/or may suffer fi'om slow performance. 

25 The discrete Fourier transform (DFT) is a popular tool for analyzing signals. 

However, before an input signal is transformed, it is quite often windowed with a 
windowing function. (It is noted that the action of capturing a finite-length sequence of 
samples of the input signal automatically implies a rectangular windowing.) The 
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frequency transform F(n) of the windowed input signal will typically exhibit multiple 
scaled and shifted versions of transform function W, i.e., the transform of the window 
function. Each sinusoidal component of the input signal expresses itself as a pair of such 
shifted versions, one version shifted up to the fi-equency ^ of the sinusoidal component, 
5 and the other shifted down to frequency -fj. The positive frequency version is referred to 
herein as a positive frequency image, and the negative frequency version is referred to 
herein as a negative frequency image. When a sinusoidal component frequency fj is small 
compared to the sample rate, the positive frequency image and the negative frequency 
image for the sinusoidal component may overlap in frequency space. Similarly, when a 

10 sinusoidal component frequency fj is close to one-half the sample rate, the positive 
frequency image and the negative frequency image for the sinusoidal component may 
overlap. Furthermore, when two sinusoidal components have frequencies that are close 
together, their positive images and negative images may overlap. 

Prior art techniques for tone estimation quite often focus on identifying the peaks 

15 in the magnitude spectrum |F(n)|. The peaks roughly determine the frequency of the 
corresponding tones. However, because of the cross-interaction of the images from other 
tones, or the negative frequency image from the same tone, the peak of a positive 
frequency image may be perturbed away from a purely scaled and frequency-shifted 
version of the template function W. Thus, parameter estimation techniques which 

20 compute parameters for a given tone based only on transform array values (i.e. DFT 
values) in the vicinity of a corresponding image peak may not produce accurate results. 
Therefore, there exists a substantial need for a system and method which can estimate 
tone parameters from the transform array for a plurality of tones with increased accuracy. 
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Summary of the Invention 

The present invention comprises various embodiments of a system and a method 
for estimating signal parameters (e.g., one or more of frequency, amplitude and/or phase) 
of a plurality of tones present in an input signal. More particularly, a system and method 
5 are described for estimating parameters for each of a plurality of tones based on a 
frequency transform F(n) of the input signal, where the input signal may be windowed 
with a window function w(n) and transformed into the frequency domain. 

In one embodiment, the method may include first receiving the input signal. For 
example, samples of the input signal may be received, where the samples include the 
10 plurality of tones. 

One or more signal parameter estimates characterizing each of the plurality of 
tones may be determined based on the received samples. In one embodiment, the one or 
more signal parameter estimates characterizing each of the plurality of tones may 
comprise a knowledge base of the signal parameter estimates, and the method may further 
15 include storing a copy of the knowledge base, as a prior knowledge base. 

Then, the method may generate a refined estimate of the one or more signal 
parameters of each respective tone of the plurality of tones, updating the one or more 
signal parameter estimates accordingly, where the refined estimates are generated in an 
iterative fashion. In other words, refined estimates of the parameters for each tone may 
20 be generated in an iterative fashion, where each successive iteration produces a more 
refined estimate of parameters for each tone. Updating the one or more signal parameter 
estimates with the refined estimate of the one or more signal parameters of the respective 
tone may include updating the knowledge base with the refined estimate of the one or 
more signal parameters of the tone, where, after updating the estimates for each tone, the 
25 knowledge base is comprised of the refined estimates for each tone. 

In one embodiment, iteratively generating the refined estimate of the signal 
parameters of each of the respective tones may include: for each respective tone of the 
plurality of tones, generating a refined estimate of the one or more signal parameters of 
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the respective tone, and updating the one or more signal parameter estimates (e.g., the 
knowledge base) with the refined estimate. The generating and updating for each 
respective tone may be performed one or more times to produce final one or more signal 
parameter estimates for each of the plurality of tones. In one embodiment, performing 
5 the generating and updating one or more times may include performing the generating 
and updating for each respective tone in an iterative manner until one or more of 1) the 
number of iterations exceeds a threshold number of iterations, and 2) the difference 
between values of the one or more signal parameters of the knowledge base and values of 
the one or more signal parameters of the prior knowledge base is less than a threshold 
10 value. 

In one embodiment, generating a refined estimate of the one or more signal 
parameters of each respective tone may include: for each respective tone of the plurality 
of tones, iteratively applying a single tone estimation method to the respective tone, 
thereby generating a refined estimate of the one or more signal parameters of the 
15 respective tone, and updating said one or more signal parameter estimates with the 
refined estimate of the one or more signal parameters of the respective tone. 

Finally, the method may include storing and/or outputting the final one or more 
signal parameter estimates for each of the plurality of tones in response to iteratively 
generating the refined estimates. 

20 

In one embodiment of the invention the method may estimate parameters for the 
plurality of tones based on a firequency transform F(n), e.g., a Fourier transform of the 
samples of the input signal. The input signal may be windowed with a window fimction 
w(n) and transformed into the fi-equency domain. In a preferred embodiment, the window 
25 fimction is a cosine-based window function. Each tone in the input signal may express 
itself in the firequency domain as an additive combination of two or more spectra, e.g., 
one centered at the tone fi-equency and the other at the negative of the tone frequency. 
These two spectra are referred to herein as the positive frequency image and the negative 
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frequency image respectively. A tone in the input signal may also be affected by spectra 
(e.g., positive or negative images) from other tones present in the signal. 

Thus, in one embodiment, the method may include generating a frequency 
transform of the samples, where the frequency transform includes frequency response 
5 data for the plurality of tones. In this embodiment, generating a refined estimate of the 
one or more signal parameters of each respective tone may include: for each respective 
tone of the plurality of tones, generating modified frequency response data for the 
respective tone which operates to remove effects of other tones of the plurality of tones 
from the frequency response data of the respective tone, using the one or more signal 

10 parameter estimates of each of the other tones, then applying a single tone estimation 
method to the modified frequency response data, thereby generating a refined estimate of 
the one or more signal parameters of the respective tone. The refined estimate of the one 
or more signal parameters of the respective tone may then be used to update the one or 
more signal parameter estimates of the respective tone. The generating, removing effects, 

15 and updating may be performed one or more times on each of the plurality of tones to 
generate final one or more signal parameter estimates for each of the plurality of tones. 

In one embodiment, removing effects of other tones from the frequency response 
data of the potential local peak of the respective tone may include: determining two or 
more sample values proximate to the potential local peak, and, for each of the other tones, 

20 calculating an effect value for each of the two or more sample values using the estimated 
signal parameters for the tone. The respective effect values from each of the two or more 
sample values may then be subtracted from the respective sample values to generate a 
corresponding two or more modified sample values, comprising the modified frequency 
response data. The modified frequency response data may then be usable in applying the 

25 single tone estimation method. In one embodiment, the two or more sample values may 
be consecutive sample values. 

Said another way, the method may transform the signal to produce frequency 
response data, then remove the effects of the other tones from each tone's frequency 
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response data to refine the parameter estimates. After each tone's parameters have been 
refined, the process may be repeated, further refining the refined estimates for each tone 
by using the refined estimates of the other tones to more accurately remove the effects of 
the other tones from the frequency response data of each tone. In other words, the refined 
5 estimates generated in each iteration may be used to remove the effects of other tones 
from the frequency response data of each respective tone, after which a further refined 
estimate may be generated for the respective tone. The method may iterate imtil the 
desired degree of refinement is achieved or some other end condition is met. 

10 In one embodiment, the frequency transform may comprise a windowed Fast 

Fourier Transform (FFT) with window function W(Py where / denotes frequency, where 
the estimated signal parameters for the tone comprise an amplitude ^4;, a frequency /y, and 
a phase ^y, and where calculating an effect value for each of the two or more sample 
values using the estimated signal parameters for the tone may comprise calculating: 

15 

dF(m) = A,e-^^W(f„ ^ f^)^ A,e-^^W(f^ -/,), 

for each sample m, wherein m denotes an index of the sample, and fm denotes the 
frequency bin of the sample, and where subtracting the respective effect value from each 
20 of the two or more sample values comprises subtracting a respective value of dF(m) from 
each of the two or more sample values. 

The plurality of tones may be processed as described one or more times to 
generate final one or more signal parameter estimates for each of the plurality of tones. 
25 Thus, the signal parameters for each of the plurality of tones may be estimated iteratively, 
updating previous parameter estimates with the refined estimates each iteration. 
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In one embodiment, the single tone estimation method applied to each tone may 
identify a frequency location proximate to an amplitude peak in the frequency transfomi, 
where the amplitude peak corresponds to the tone. Two or more frequency bins may then 
be selected proximate to the identified frequency location in the frequency transform, 
5 The method preferably selects one or more frequency bins located on either side of the 
frequency location of the amplitude peak. In other words, the method may select one or 
more bins located on one side of the amplitude peak and one or more bins located on the 
other side of the amplitude peak. 

The method may then determine a tone frequency value that minimizes a 
10 difference between two or more expressions, e.g., at least a first expression and a second 
expression. Stated another way, the method may select a tone frequency value that makes 
a plurality of expressions most nearly equal. 

Each of the plurality of expressions may comprise a sum of one or more 
numerator terms divided by a sum of one or more denominator terms, wherein: 
15 1) each of the plurality of expressions includes a tone frequency variable; 

2) each numerator term and each denominator term corresponds to one of the 
frequency bins; and 

3) a ratio of each numerator term and its corresponding denominator term 
represent a complex amplitude of the tone at a respective bin. 

20 Thus, for example, the plurality of expressions may comprise a first expression 

and a second expression having the above characteristics. 

In each expression, the tone frequency variable may represent a correct tone 
frequency value of the tone. The first expression may be approximately equivalent to the 
second expression when the correct tone frequency value is used for the tone frequency 
25 variable in the first and second complex expressions. 

The step of determining a tone frequency value may comprise computing a 
plurality of differences between the first expression and the second expression for 
different respective tone frequency values of the tone frequency variable, and then 
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selecting the tone frequency value that produces a smallest difference. The method may 
compute the pluraUty of differences by performing a Newton-Rhapson root finding 
method. 

The expressions may be real expressions or complex expressions. Where the 
expressions are complex expressions, the method may involve minimizing a difference 
between an amplitude of the first complex expression and an amplitude of the second 
complex expression. 

In a preferred embodiment, the first expression and the second expression have 
the form: 

k+M . X 

'Z{F(n)W\f„ -f,)-F\n)W{f„ 
'f\w(J„-f,f-\W{f„^f,f) 

wherein: 

F(n) is the wth value of the single sided scaled FFT spectrum; and 
W represents a window function, wherein the window fimction is shifted by a 
value of the tone frequency variable fj. 

When the frequency transform of the samples comprises generating a power 
spectrum of the samples, the first expression and the second expression have the form: 

-/•)!) 

_M 

Jfc+1 

wherein: 



tXnnW\L -L)-F\n)W{f„ +/.)) 

ifc+A/-l / \ 

Jfc+1 



ifc+1 

Jfc+A/-1 

Z \w{f„-fd\ 

ifc+1 
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F(n) is the nth value of the single sided scaled FFT spectrum; and 
W represents a window function, wherein the window function is shifted by a 
value of the tone frequency variable fj. 



correct frequency of the tone. In other words, the determined tone frequency value 
represents an approximation of the correct tone frequency value. The method described 
herein operates to find a more exact approximation of the frequency of the tone than prior 
art techniques, even in the presence of noise or interference from other tones, or the 
10 negative image of the respective tone being found. 

After the tone frequency value has been determined, the method may comprise 
storing the determined tone frequency value in a memory, or outputting the determined 
tone frequency value, e.g., on a display or to other systems for fiirther analysis or 
processing. 

15 The method may also then compute one or more of the amplitude and phase of the 

tone using the determined tone frequency value. For example, once thefj value is known, 
the amplitude and phase can be computed as: 



5 



The determined tone frequency value may comprise an approximation of the 



A.e'^^ =C = 




20 



When the frequency transform of the samples comprises generating a power 
spectrum of the samples, the amplitude can be computed as: 
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4= 



^(\Fin)\x\W(f„-f,)\) 



z 



Ei^(/--/)r 



As mentioned above, the above steps may be iteratively perfomied for each of the 
plurality of tones to determine at least one signal parameter for each of the plurality of 
5 tones in the input signal. The above method may be performed iteratively where, for 
each respective tone identified, the other tones' effects may be subtracted from the input 
signal, and the above single tone estimation process performed on the signal The process 
may then be repeated for the next tone, and so on, generating refined estimates for the 
parameters of each tone. 

10 As also described above, once all of the tones have been refined by this 

subtraction, the entire process may be repeated (iterated), using the refined estimates to 
subtract the effects of the other tones from the signal, then performing single tone 
estimation on the signal. Each process iteration may generate a fiirther refined set of 
parameter estimations for the tones. The process iteration may be terminated upon 

15 reaching a threshold number of iterations, or when the difference between successive 
estimates of the tone parameters is less than a threshold value, i.e., when the desired 
degree of refinement has been achieved, resulting in a final set of estimates for the tone 
parameters. The final estimates of the tone parameters may then be stored, or output to a 
display or external system. 

20 It is noted than in other embodiments, altemative single tone estimation methods 

may be used in place of the above single tone estimation method, such as, for example, 
the method described in U.S. Patent Application Serial No. 09/753,164, titled "System 
and Method for Estimating Tones in an Input Signal", which is incorporated by reference 
below. 

25 



Atty.Dkt.No.: 5150-62101 



Page 10 Meyertons, Hood, Kivlin, Kowert & Goetzel PC 



Thus the method of the present invention operates to determine more accurate 
approximations of signal parameters of multiple tones present in an input signal, and 
operates more efficiently than prior art techniques. 
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Brief Description of the Drawings 

A better understanding of the present invention can be obtained when the 
following detailed description of the preferred embodiment is considered in conjunction 
with the following drawings, in which: 
5 Figure lA illustrates a system configuration 100 for determining the signal 

parameters associated with one or more sinusoidal tones comprised within an input 
signal; 

Figure IB illustrates one embodiment for tone detection system 120; 
Figure 2 A illustrates one embodiment of tone detection system 120 comprising a 
10 computer-based measurement or data acquisition system, where signals generated by a 
signal reception device (SRD) are presented to computer 102 through signal conditioning 
system 108 and data acquisition (DAQ) device 104; 

Figure 2B illustrates a second embodiment of tone detection system 120 
comprising a computer-based measurement system, where signals generated by a signal 
15 reception device (SRD) are presented to computer system 102 through data acquisition 
(DAQ) device 104; 

Figure 3 illustrates peak magnitudes of transform array F(k) for a typical 
windowed input signal comprising multiple sinusoidal tones. 

Figure 4 flowcharts a method for estimating parameters for multiple tones in a 
20 signal, according to one embodiment; 

Figure 5 flowcharts a more detailed embodiment of the method of Figure 4; and 

Figures 6A and 6B, 7A and 7B, 8A and 8B, and 9A and 9B illustrate front panels 
and block diagrams of a portion of a graphical program that implements one embodiment 
of the method of Figure 5. 
25 Figure 10 presents a flowchart for one embodiment of a single tone detection 

method according to the present invention; 

Figure 1 1 illustrates a windowing operation being performed on an input signal to 
generated a windowed input signal; 



1 
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Figure 12 illustrates the magnitude of transform array F(k) for a typical windowed 
input signal comprising a single sinusoidal tone; 

Figure 13 illustrates a blowup of a generic magnitude peak from the magnitude 
spectrum of Figiwe 12; 

5 Figures 14A and 14B, 15A and 15B, 16A and 16B, 17A and 17B, and 18A and 

18B illustrate front panels and block diagrams of a portion of a graphical program that 
implements one embodiment of the invention. 

While the invention is susceptible to various modifications and altemative forms, 
10 specific embodiments thereof are shown by way of example in the drawings and will 
herein be described in detail. It should be understood, however, that the drawings and 
detailed description thereto are not intended to limit the invention to the particular form 
disclosed, but on the contrary, the intention is to cover all modifications, equivalents and 
altematives falling within the spirit and scope of the present invention as defined by the 
1 5 appended claims. 



Atty. Dkt. No.: 5150-62101 



Page 13 



Meyertons, Hood, Kivlin, Kowert & Goetzel PC 



Detailed Description of the Embodiments 

Incorporation by Reference 

U.S. Patent Application Serial No. 09/753,164 titled "System and Method for 
5 Estimating Tones in an Input Signal", filed December 27, 2000, and whose inventor is 
Alain Mariot, is hereby incorporated by reference in its entirety as though fully and 
completely set forth herein. 

U.S. Patent AppUcation Serial No. 09/880,167 titled "System and Method for 
Estimating One or More Tones in an Input Signal", filed June 12, 2001, and whose 
10 inventors are Yong Rao and Michael Cema, is hereby incorporated by reference in its 
entirety as though fiiUy and completely set forth herein. 

Figure lA 

Figure 1 A illustrates a system configuration 100 for performing signal processing 
15 on a signal comprising a plurality of tones. System configuration 100 may comprise a 
signal reception device (SRD) and a tone detection system 120. 

The SRD may be coupled to receive a signal firom a device, unit under test (UUT) 
or a transmission medium 110, or any other system capable of transmitting a signal that 
may contain tones. The term "transmission medivun" is used herein to refer generally to a 
20 ' device, unit under test (UUT) or a transmission medium 110 that may generate a signal 
including two or more tones. As used herein, the term "tone" includes a signal at a 
firequency, e.g., at a primary or single fi-equency, which may be contained within another 
signal or may itself be the signal. For example, the signal may comprise two or more 
tones and a small amount of noise. As an example, a tone may comprise a sinusoidal 
25 signal having arbitrary amplitude, fi-equency and phase. 

As shown in Figure lA, the SRD may be coupled to a transmission medium 110. 
Transmission medium 110 may represent any of a variety of transmission media such as 
the atmosphere, fi-ee space, an optical fiber or fiber bundle, a communication bus (e.g. a 
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network bus), a body of water or any other fluid, the earth, etc. In one embodiment, the 
transmission medium 110 is the atmosphere, and the SRD comprises an antenna and a 
radio receiver. In a second embodiment, the transmission medium 110 is a network bus 
connecting two or more computers, and the SRD is a network interface card/board. In a 
5 third embodiment, the transmission medium 110 is an optical fiber, and the SRD 
comprises an optical sensor. As noted above, element 110 may be any of various devices 
or mediums for generating or transmitting a signal 

The SRD receives an input signal from the transmission medium or device 110 
and may convert the input signal into a form suitable for presentation to tone detection 
10 system 120. The input signal may be electrical or non-electrical in nature. The SRD may 
include analog-to-digital conversion hardware to digitize the input signal. Altematively, 
analog-to-digital conversion hardware may be comprised within tone detection system 
120. 

In one embodiment, the SRD may comprise a measurement device such as a 
15 microphone, an accelerometer, a spatial displacement sensor, a strain gauge, a pressure 
sensor, a temperature sensor (e.g., a thermocouple), a radiation sensor, an optical sensor, 
etc, or any combination thereof. In another embodiment, the SRD may represent an array 
of transducers or measurement devices of one or more types. The SRD may thus be any 
of various transducers or sensors for receiving a signal. 
20 Tone detection system 120 may couple to the SRD to receive the input signal. In 

an alternate embodiment, the tone detection system 120 may be operable to receive the 
input signal from a memory, wherein the input signal was previously acquired and stored 
in the memory. 

Tone detection system 120 may be configured for detecting the frequency, 
25 amplitude and/or phase of a plurality of tones in the input signal. Tone detection system 
120 may comprise a processor or central processing unit (CPU) 140, memory 146, user 
input device(s) UID and a display device DD as shown in Figure IB. CPU 140 may be 
realized by any of a variety of computational devices such as a general pvupose processor. 
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a digital signal processor, a parallel processor, dedicated digital and/or analog circuitry, 
programmable gate array logic (e.g., an FPGA), etc., or any combination thereof. 
Memory 146 may comprise any of a variety of memory devices such as random access 
memory (RAM) and/or read-only memory (ROM), as described further below. Tone 
5 detection system 120 may also include specialized data acquisition and/or signal 
conditioning hardware, interface hardware, etc., or any combination thereof 

Tone detection system 120 may comprise any of various devices, such as a 
programmable computer system, a computer-based system such as a VXI-based system, a 
PXI-based system, a GPIB-based system, a computer-based data acquisition system, or a 

10 dedicated test instrument, such as a dynamic signal analyzer, an oscilloscope or any other 
signal acquisition and/or analysis device. 

Tone detection system 120 may operate on samples of the input signal X 
generated by the SRD and thus may identify the frequency, phase and/or amplitude of 
multiple tones in the input signal. The frequency, phase and/or amplitude of the tones 

15 may be presented to a user through the display device DD or some other output device, 
may be further processed by a software program executing in the tone detection system 
120, may be provided to another system for processing, and/or may be stored to memory 
for future use. 

User input device(s) UID may comprise a keyboard, a pointing device such as a 
20 mouse or trackball, a touch pad (such as those used in modem laptop computers for 
cursor control), a touch sensitive display screen, etc., or other input devices. In one 
embodiment, user input device(s) UID may include use of a graphical user interface or 
panel configured with various control icons such as buttons, knobs, sliders, switches, 
indicators, etc., or any combination thereof A user provides input to tone detection 
25 system 120 through user input device(s). Tone detection system 120 may manage a 
graphical user interface through display device DD and user input device(s) UID. 

Figures 2A and 2B 
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Figure 2A and 2B illustrate exemplary embodiments of tone detection system 
120. As shown, tone detection system 120 may comprise a computer 102, a data 
acquisition PAQ) device 104 coupled to the computer 102, and optionally a signal 
conditioning system 108 coupled to the DAQ device 104. The SRD may comprise 
5 transducers, sensors, and/or receiving devices ttiat couple to DAQ device 104 through the 
signal conditioning circuitry 108. 

As shown, the SRD is configured and/or coupled to acquire signals from the 
transmission medium 110. The input signals acquired by the SRD may be optionally 
conditioned by the signal conditioning system 108 as shown in Figure 2 A. The conditioned 

10 input signals may then be provided to DAQ device 104 as shown. Signal conditioning 
system 108 may connect to DAQ device 104 via one or more cables. 

In one embodiment, signal conditioning system 108 may comprise an extemal 
chassis 122 housing one or more signal conditioning modules 124 and optionally terminal 
blocks 126. Signal conditioning system 108 may be used to perform signal conditioning on 

15 the signals generated by the SRD. As used herein, the term "signal conditioning" may 
include one or more of amplifying, linearizing, limiting, isolating, filtering, switching 
and/or multiplexing signals), among other signal processing functions. Signal 
conditioning system 108 may advantageously reduce the level of noise in the signals 
transmitted to DAQ device 104. DAQ device 104 may receive conditioned signals from 

20 signal conditioning system 108 as shown in Figure 2 A. Alternatively, DAQ device 104 
may directly receive the input signal from the SRD as shown in Figure 2B. DAQ device 
104 may operate to perform analog to digital (A/D) conversion and provide the resultant 
digital signals to computer 102 for processing. 

Computer system 102 may include various standard components, including a 

25 processor or central processing unit (CPU) 140, system memory 146, non- volatile 
memory, one or more buses, and a power supply. DAQ device 104 may be a specialized 
system for acquiring digital and/or analog signals from extemal devices. Thus, DAQ 
device 104 may include analog to digital (A/D) conversion circuitry and/or digital to 
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analog (D/A) conversion circuitry. Examples of the DAQ device 104 include "E series" 
DAQ boards from National Instruments Corporation. DAQ device 104 may also 
comprise a computer-based instrument board, such as an oscilloscope, a digital multimeter 
(DMM), a dynamic signal analyzer, an arbitrary waveform generator, etc. 
5 In one embodiment, computer 102 may comprise input/output (I/O) slots into 

which DAQ device 104 may be coupled. In another embodiment, computer 102 may 
comprise a VXI (VME Extensions for Instrumentation) chassis and bus, a GPIB (General 
Purpose Interface Bus) interface card, a serial port or parallel port by which DAQ device 
104 may be coupled to the computer 102. 

10 Tone detection system 120, e.g., computer system 102, preferably includes at least 

one memory medium on which computer programs according to the present invention may 
be stored. The term "memory medium" is intended to include various types of memory or 
storage, including an installation 'medium, e.g., a CD-ROM, or floppy disks 104; a 
computer system memory or random access memory such as DRAM, SRAM, EDO RAM, 

15 Rambus RAM, EPROM, EEPROM etc.; or a non- volatile memory such as a magnetic 
media, e.g., a hard drive, or optical storage. The memory medium may comprise other 
types of memory as well, or combinations thereof. In addition, the memory medium may 
be located in a first computer in which the programs are executed, or may be located in a 
second different computer which connects to the first computer over a network. In the latter 

20 instance, the second computer may provide the program instructions to the first computer 
for execution. Also, the computer system 102 may take various forms, including a personal 
computer system, mainfii^e computer system, workstation, network appliance, Internet 
appliance, personal digital assistant (PDA), television system, dedicated test or 
measurement instrument or other device. In general, the term "computer system" can be 

25 broadly defined to encompass any system having a processor which executes instructions 
from a memory medium. 

The memory medium preferably stores a software program according to one 
embodiment of the present invention for detecting multiple tones in the input signal. 
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More particularly, the software program may be operable to analyze the input signal to 
determine the frequency, phase and amplitude of multiple tones in the input signal. 

The software program may be implemented in any of various ways, including 
procedure-based techniques, component-based techniques, object-oriented techniques, or 
5 neural net based learning techniques, among others. For example, the software program 
may be implemented using ActiveX controls, C-H- objects, Java objects, Microsoft 
Foundation Classes (MFC), or other technologies or methodologies, as desired. A 
processor, such as the host CPU, executing code and data from the memory medium, or a 
programmable device configured according to a net list, may comprise embodiments of a 

10 means for determining the frequency, phase and amplitude of the tones embedded in the 
input signal according to the methods described below. 

As noted above, the tone detection system 120 may also or instead include 
reconfigurable logic or programmable logic, such as an FPGA, that implements the 
methods described herein. 

15 Various embodiments fiirther include receiving, storing, and/or transmitting 

instructions and/or data implemented according to the present invention upon a carrier 
medium. Suitable carrier media include a memory medium as described above, as well as 
signals such as electrical, electromagnetic, or digital signals, conveyed via a 
communication medium such as networks and/or a wireless link. 

20 

Multiple Tone Estimation 

As mentioned above, in many appUcations, the received signal includes a plurality 
of tones. Various embodiments of a method for estimating the signal parameters for each 
of a plurality of tones in a received signal are described below with reference to Figures 
25 3-9B. In particular, various embodiments of a method are described which may be used 
to identify and estimate signal parameters of a plurality of tones in a signal, compensating 
for the effects of both the negative image of each tone, and the positive and negative 
images of the other tones present in the signal. In estimating the signal parameters for the 
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plurality of tones, the method preferably utiHzes a single tone estimation method. 
Exemplary embodiments of such a single tone estimation method are described below 
with reference to Figures 10-1 8B. 

5 Figure 3 - Frequency Response Data With Multiple Tones 

Figure 3 is a plot of the magnitude of a transform F(n) corresponding to a typical 
windowed input signal x(n), illustrating frequency response data generated by the 
frequency transform F(n) of a signal containing multiple tones, where the horizontal k 
axis denotes frequency. As Figure 3 shows, each tone is represented by a pair of peaks 

10 symmetrically located about the origin — one at the positive frequency, and another at the 
negative frequency. The positive peak for each tone may be detected and characterized 
using the transform values in the neighborhood of, or proximate to, the peak. 

It is noted that in defining a peak, e.g., peak 301, the maximum of magnitude 
fimction |F(n)| considered as a ftmction of continuous frequency typically does not occur 

15 at the integer value k^^, although it should occur somewhere in the interval between k^^^x 
and the next highest amplitude sample adjacent to k^^. A comparison of amplitude at 
k^-1 and k^+1 may be made to determine whether the second largest magnitude occurs 
at kn^-1 or k^+1. Let kz denote the location of this second largest magnitude. Figure 3 
illustrates a blowup of the positive frequency magnitude peak 301 in the neighborhood of 

20 bin index value k^^, and a second positive frequency magnitude peak 303 in the 
neighborhood of bin index value k'^. In the example of Figure 3, peak 301, the 
frequency bins k^^-l and kj^-^l may be used as the bins proximate to the frequency 
producing the positive maximum amplitude. Although in this example, three data points 
are used to define a peak, it is noted that in other embodiments two points may be used to 

25 define each peak, or altematively, four or more points may be used. 

Figure 4 - Flowchart of a Method for Multi-Tone Estimation 
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Figure 4 flowcharts one embodiment of a method for detecting multiple tones in a 
signal. It is noted that in other embodiments, some of the steps may occur in a different 
order than shown, or may be omitted. Additional steps may also be performed as desired. 
In a preferred embodiment, the method may be implemented by one or more tone 
5 detection software programs executed by a processor, e.g., in a computer system. The 
programs may be stored in a memory medium coupled to the processor. In another 
embodiment, the method may be implemented in programmable hardware, such as an 
FPGA. In further embodiments, the method may be implemented by a combination of 
software and hardware. 

10 As Figure 4 shows, in 1302 an input signal may be received which contains two 

or more tones. In other words, samples of the signal may be received, such as through an 
input of the computer system, where the signal samples include a plurality of tones. 

In 1304, a frequency transform of the signal samples may be generated, resulting 
in frequency response data for the plurality of tones. In one embodiment, the frequency 

15 transform is a windowed Fast Fourier Transform, i.e., a windowed Discrete Fourier 
Transform, as described below with reference to Figures 10 and 11. As is known in the 
art, the frequency response data includes amplitudes corresponding to each frequency 
value in the signal spectrum. 

In 1306, a plurality of potential local peaks may be determined based on the 

20 frequency response data, where each potential local peak corresponds to a tone in the 
signal. Said another way, a plurality of potential local maxima (positive peaks) may be 
determined based on the frequency transform. In one embodiment, the peaks or maxima 
may be determined by locating sample amplitudes exceeding a predetermined threshold. 
It should be noted that the term "potential" is used because each potential peak is defined 

25 by two or more discrete samples or data points, and the actual peak amplitude generally 
occurs between two of the samples, as indicated in Figure 5. In other words, it is unlikely 
that one of the discrete samples lies precisely at the peak or local maxima, thus, the 
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sample values defining each peak indicate a potential peak, which may be refined as 
described below. 

In 1308, one or more signal parameter estimates characterizing each of the 
plurality of tones may be determined based on the plurality of potential local peaks. 
5 These estimates may be considered initial estimates which provide a rough starting point 
for the estimation process described below. In one embodiment, the one or more signal 
parameter estimates may include estimates for one or more of firequency /, amplitude Ay 
and phase (|), for each tone. In a preferred embodiment, the estimate of firequency may be 
the fi-equency value of the determined potential local peak, i.e., the firequency value of the 

10 maximum sample amplitude for that local peak. The initial estimate of amplitude may be 
zero. Similarly, the initial estimate of phase may be zero. Thus, for each potential local 
peak, i.e., for each tone, initial estimates may be determined which characterize the tone 
associated with each potential local peak in the fi-equency response data. 

In 1310, the one or more signal parameter estimates characterizing each of the 

15 plurality of tones may be stored, e.g., in the memory medium of the computer system. In 
one embodiment, storing the one or more signal parameter estimates for each of the 
plurality of tones may include sorting the signal parameter estimates for the plurality of 
tones by the amplitude of the potential local peaks. In other words, the tone parameters 
may be sorted by descending peak amplitude. 

20 Steps 1314 and 1316 may then be performed for each estimated tone, i.e., for each 

tone whose initial parameter estimates have been determined and stored, as shown in 
Figure 4 and described below. 

For each respective tone, in 1314, the effects of the other tones of the plurality of 
tones may be removed fi"om the fi-equency response data of the tone. Said another way, 

25 the influence of each of the other tones may be subtracted fi:'om the samples defining the 
respective tone's potential local peak. In one embodiment, removing the effects of other 
tones fi-om the firequency response data of the potential local peak of the tone includes 
determining two or more sample values proximate to the potential local peak, then, for 
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each of the other tones, calculating an effect value for each of the two or more sample 
values using the estimated signal parameters for the tone. The respective effect values 
may then be subtracted from each of the two or more sample values, thereby generating a 
corresponding two or more modified sample values. Thus, the original sample data 
5 defining the potential local peak may be modified to produce new sample data with the 
effects of the other tones removed. 

As mentioned above, in one embodiment, the estimated signal parameters for each 
tone include an amplitude a frequency //, and a phase (pj, and the frequency transform 
is a windowed Fast Fourier Transform (FFT) with window fimction W(f), where / denotes 
10 frequency. In this embodiment, the effect value of each of the tones may be calculated 
for each of the two or more sample values by the following expression: 



15 where m denotes an index of the sample, and^ denotes the frequency value or bin 

of the sample. As noted above, the amplitude Aj^ frequency //, and phase ^y, are the 
estimated parameters used to characterize each of the other tones. Thus, subtracting the 
respective effect value of each tone from each of the two or more sample values includes 
calculating and subtracting a respective value of dF(m) from each of the two or more 

20 sample values, per tone to be removed. Thus, the method may iterate through each of the 
other tones, calculating and removing each tone's effect from each of the samples. It is 
noted that in the first pass through the process, the initial estimated amplitude Ai^ and 
phase (f>i, may be zero, and so the calculated effect value dF(m) for each of the other tones 
may be zero for the first respective tone. In an embodiment where the stored tone 

25 parameters are sorted (e.g., base on peak amplitude), step 1314 may process the tones in 
the sorted order. This may increase the efficiency of the estimate refinement process. 
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After the removal of effects of the other tones in 1314, a single tone estimation 
method may be applied to the modified frequency response data to generate a refined 
estimate of signal parameters for the tone, as indicated in 1316. The refined estimate of 
the signal parameters of the tone may replace the original estimate of the signal 
5 parameters, and thus may be used in subsequent calculations of the effect value of the 
tone, e.g., when calculating the effect of the tone on each of the other tones. It is noted 
that the single tone estimation method of 1316 preferably comprises the single tone 
estimation method described below with reference to Figures 10-13, or variations thereof 
As indicated in 1320, if there are more tones to be estimated, then the method may 

10 continue with 1314 above, performing steps 1314 and 1316 for each tone, thereby 
generating a refined estimate of the signal parameters for each of the tones in the signal. 

In one embodiment, if in 1320, it is determined that there are no more tones to be 
processed, then the method may repeat steps 1310 through 1320 one or more times, 
storing the refined estimates of the signal parameters of the tones, using the refined 

15 estimates to remove the effects of the other tones from each tone, and applying the single 
tone estimation method to each tone to generate fiuther refined estimates. In other words, 
the refinement process may be repeated one or more times, where the removal of effects 
performed in 1314 may be more accurate with each succeeding iteration. Further details 
of such an iterative approach to estimation refinement are presented below with reference 

20 to Figure 5. 

Figure 5 - Detailed Flowchart of a Method for Multi-Tone Estimation 

Figure 5 is a detailed flowchart of another embodiment of a method for detecting 
multiple tones in a signal. As mentioned above, in various embodiments, some of the 
25 steps may occur in a different order than shown, or may be omitted. Additional steps 
may also be performed as desired. As may be seen, several of the steps described below 
are substantially the same as corresponding steps of the method of Figure 4, and 
descriptions thereof may be abbreviated. 
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As Figure 5 shows, in 1302 an input signal may be received which contains two 
or more tones. In other words, samples of the signal may be received, where the signal 
samples include a plurality of tones. 

In 1304, a frequency transform of the signal samples may be generated, resulting 
5 in frequency response data for the plurality of tones. The frequency transform is 
preferably a windowed Fast Fourier Transform, i.e., a windowed Discrete Fourier 
Transform, as described below with reference to Figures 10 and 11. 

In 1306, a plurality of potential local peaks may be determined based on the 
frequency response data, where each potential local peak corresponds to a tone in the 
10 signal. 

Then, in 1309, a knowledge base may be built based on the pluraUty of potential 
local peaks. In one embodiment, the knowledge base may include one or more signal 
parameter estimates characterizing each of the plurality of tones. As mentioned above in 
1308, these estimates may be considered initial estimates which provide a rough starting 

15 point for the estimation process described below, and may include estimates for one or 
more of frequency /, amplitude A, and phase (|), for each tone. As also mentioned above, 
in a preferred embodiment, the frequency estimate may be the frequency value of the 
determined potential local peak, i.e., the frequency value of the maximum sample 
amplitude for that local peak, and the initial estimates of amplitude and/or phase may be 

20 zero. Thus, a knowledge base may be generated which contains initial estimates for one 
or more signal parameters for each tone. 

In one embodiment, the tone parameter estimates in the knowledge base may be 
sorted. For example, the tone parameters may be sorted by amplitude of the 
corresponding potential local peaks. Subsequent processing of the tones may then be 

25 performed in the sorted order. 

In 1311, a copy of the knowledge base built in 1309 may be stored as a prior 
knowledge base, e.g., in the memory medium of the computer system. In other words, a 
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copy of the knowledge base may be saved for later use by the method, where the stored 
copy comprises a "prior version" of the knowledge base. 

Then, as Figure 5 shows, steps 1315, 1316, and 1318 may be performed for each 
estimated tone, i.e., for each tone whose parameter estimates have been determined and 
5 stored in the knowledge base. Thus, for each respective tone represented in the 
knowledge base, in 1315, the effects of the other tones of the plurality of tones may be 
removed from the frequency response data of the tone, similar to step 1314 of Figure 4 
above. 

After the removal of effects of the other tones in 1315, a single tone estimation 
10 method may be applied to the modified frequency response data to generate a refined 
estimate of signal parameters for the tone, as previously indicated in 1316. As noted 
above, the single tone estimation method of 1316 preferably comprises the single tone 
estimation method described below with reference to Figures 10-13, or variations thereof 
Then, in 1318, the knowledge base may be updated with the refined estimate 
15 generated in 1316. It is noted that this update of the knowledge base preferably does not 
modify the stored "prior" knowledge base from 1311. Maintaining a prior version of the 
knowledge base allows comparison of the two versions to assess the need for fiirther 
refinement of the estimates, as described below with respect to the determination of end 
conditions. 

20 As indicated in 1320, if there are more tones to be estimated, then the method may 

continue with 1315 above, performing steps 1315, 1316, and 1318 for each tone, thereby 
generating an updated knowledge base containing the refined estimates of the signal 
parameters for each of the tones in the signal. 

If in 1320, it is determined that there are no more tones to be processed, then in 

25 1322, a determination may be made as to whether an end condition is met. If the end 
condition is met, the method may terminate, as shown. It is noted that this conditional 
controls the iteration of steps 1311, 1315, 1316, 1318, and 1320, described above. 
Examples of the end condition are presented below. Thus, the method may continue with 
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step 1311, storing the updated knowledge base of refined estimates as a new "prior" 
knowledge base, i.e., replacing the old prior knowledge base. The method may then 
repeat steps 1311-1320, until the end condition of 1322 is met, using each successive 
updated knowledge base to remove the effects of the other tones from each tone, and 
5 applying the single tone estimation method to each tone to generate further refined 
estimates. In other words, the refinement process of steps 1315-1320 may be repeated 
one or more times. Because the updated knowledge base is used to remove the effects of 
the other tones from each tone, and the knowledge base is updated every iteration, the 
removal of effects performed in 1315 may be more accurate with each succeeding 
10 iteration. 

In one embodiment, the end condition of 1322 may comprise the number of 
iterations meeting or exceeding a threshold number, e.g., four. In another embodiment, in 
1322, the difference between values of the signal parameters of the updated knowledge 

15 base and values of the signal parameters of the prior knowledge base may be calculated, 
and if the difference is less than a threshold value, the iteration may terminate. In yet 
another embodiment, either of both conditions may terminate the iteration. Thus, if the 
refinement process of steps 1315-1318 do not change or improve the parameter values of 
the knowledge base more than the threshold value, the method may terminate. In one 

20 embodiment, the latest version of the knowledge base, i.e., the latest version of the 
refined signal parameters of the tones, may be output, where the refined signal parameters 
comprise final estimates of the signal parameters of the tones. 

Example Graphical Program for Multi-Tone Detection 
25 Figures 6A -9B comprise screen shots of a portion of a graphical program which 

implements one embodiment of the invention. Figures 6 A -9B comprise screen shots of a 
portion of a graphical program or Virtual Instrument (VI) written in the Lab VIEW 
graphical programming language. The graphical program is hierarchical in nature, where 
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Figures 6A and 6B comprise the top level front panel and block diagram, and figures 
7A/B, 8A/B, and 9A/B are various sub programs or sub- Vis of the program. Various 
other sub programs or sub-Vis of this graphical program that are not necessary to an 
understanding of this embodiment of the invention are not included. 
5 Figures 6A and 6B illustrate a top level front panel and block diagram, 

respectively, of the graphical program titled "multi-tone estimation example". As may be 
seen in Figure 6A, this panel provides a top level interface for the user to configure the 
multi-tone estimation process, as well as to view the results of the process. As Figure 6B 
shows, this graphical program block provides top level control and execution of the 

10 multi-tone estimation process. This graphical program includes a block titled "spectrum 
peak tone search and estimation (FFT spectrum)''. 

Figures 7A and 7B illustrate a front panel and block diagram, respectively, of the 
graphical program block titled "spectrum peak tone search and estimation (FFT 
spectrum)" contained in the diagram of Figure 6B. This graphical program provides the 

15 primary fimctionality of the invention, including blocks titled "spectrum compensation 
from other tones" and "SML tone estimation (FFT input)". 

Figures 8A and 8B illustrate a front panel and block diagram, respectively, of the 
graphical program block titled "spectrum compensation from other tones" contained in 
the diagram of Figure 7B. As the title indicates, this graphical program calculates and 

20 removes the effects of the other tones from a given tone. 

Figures 9 A and 9B illustrate a front panel and block diagram, respectively, of the 
graphical program block titled "SML tone estimation (FFT input)" contained in the 
diagram of Figure 7B. This graphical program performs the single tone estimation 
method on the modified or compensated frequency response data of a tone, and includes 

25 two blocks labeled "G(f)" and "G'(f)"» whose front panels and block diagrams are shown 
in Figures 17A and 17B, and Figures ISA and 18B, respectively, and described below. 

Figure 10 - Determining Signal Parameters of a Tone 
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Figure 10 illustrates one embodiment of a method for performing single tone 
estimation of one or more signal parameters, e.g., one or more of the frequency, 
amplitude and/or phase, of a tone present in the input signal. In particular, the method 
illustrated in Figure 10 is preferably used to perform the single tone estimation referred to 
5 in 1316 of Figures 4 and 5 above. In one aspect, the method may be referred to as an 
aliasing compensation method for compensating for the effects of the negative image of a 
tone on the positive image of the tone. In another aspect, the method may be used to 
more correctly identify signal parameters of a tone in a signal, compensating for the 
effects of the negative image of the tone present in the signal. 

10 The method of Figure 10 may be implemented by execution of a computer 

program stored on the memory medium as described above. The method may also be 
implemented in programmable logic, such as an FPGA. The method may also be 
implemented by a combination of hardware and software, as desired. It is noted that 
various steps may be performed concurrently or in a different order than that shown, 

15 and/or some steps may possibly be omitted, as desired. For example, when the method of 
Figure 10 is used for the single tone estimation in step 1316 of the methods of Figures 4 
and 5, the steps 202 and 204 below may already have been performed, and so may be 
omitted. 

In step 202 the method may receive samples x(n) of the input signal, wherein the 
20 input signal includes the tone (or multiple tones). The samples x(n) may be provided by 
the SRD or may be received from a memory medium, e.g., having been previously 
recorded/captured from the SRD. Altematively, the input signal samples may be 
simulated samples generated by a simulator (e.g. a CPU executing simulation code). 
Various embodiments of the invention contemplate a wide variety of possible sources for 
25 the input signal samples x(n). 

The input signal may comprise a single sinusoidal tone in the presence of noise. 
Thus, the input signal and/or the tone may be modeled by the expression: 
x(n) = i4. cos(2;^.n + <p-) 
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Where Ai is the amplitude of the tone,/ is the normaUzed frequency of the tone, n 
is from 0 to AT-i, and ^ is the initial phase of the tone. 

In step 204 the CPU 140 may multiply the input samples by a known window 
5 function w(n) to generate a windowed input signal y(n) = w(n) * x(n) as suggested by 
Figure 11, Alternatively, the CPU 140 may receive, e.g., in step 202, a discrete set of 
samples, which effectively form a windowed version of the signal. The window fimction 
w(n) may have any of a variety of forms. For example, the window function may be a 
rectangular window, a triangular window, a raised cosine window, a Hanning window, 
10 etc. The window function may, for example, be any type of cosine window, such as a 
Hanning, Hamming, Blackman Harris, and Flat Top window, among others. 

In step 204, the method may also perform a frequency transform, such as the 
discrete Fourier transform (DFT), e.g., an FFT (Fast Fourier Transform), on the 
windowed input signal y(n) to generate a transform array F(n), where n is a frequency bin 
15 index which may range from 0 to N-1, or any interval of length N, where N is a positive 
integer. In one embodiment or one example, the transform array F(n) may be modeled by 
the transform of the sinusoidal tone, i.e., 

F(n) = A,e^^^ W(f„ + ) + A^e^'^^ W(f„ - ) 
where W(f) represents the Fourier transform of the window w(n). It is noted that the 
20 relationship between frequency f and frequency bin number k is given by 

f=fs*(k/N), 

where fg is the sample rate. The magnitude of the window transform W(f) typically has 
even symmetry and attains a maximum at f = 0. Thus, the function W(f-fi) attains a 
maximum magnitude at frequency f = f;, and the function W(fH-fi) attains a maximum 
25 magnitude at frequency f = -f^. The first term in the expression above, i.e., 

P(f) = (A/2)exp(j^)W(f-fi), 
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is referred to herein as the "positive-frequency image" since its center frequency occurs at 
the positive frequency fj. The second term in the expression above, i.e., 

N(f) = (A/2)exp(-j^)W(f-ff;) 
is referred to herein as the "negative-frequency image" since its center frequency occurs 
5 at the negative frequency -fj. Thus, the transform array F(n) may include a positive- 
frequency image and negative-frequency image which combine additively (in the sense of 
complex addition). The input signal may also include noise, other spurious tones, or 
other valid tones. 

If tone frequency fj stays away from zero or fs/2, and/or the sample size N is 
10 sufficiently large, the overlap between the positive and negative frequency images may be 
small, and thus, their individual identities may be apparent in the transform array F(n). 
The magnitude fimction |F(n)| will thus exhibit two peaks which correspond to the 
positive and negative frequency images. The frequency locations of one of these peaks 
(i.e., the peak that occurs in the range of positive frequencies) may be used as an estimate 
1 5 for the tone frequency f^. 

Conversely, if the tone frequency is close to zero or fs/2, and/or, the sample size N 
is sufficiently small, the positive-frequency image and negative frequency image may 
overlap significantly. Thus, their individual identities may not be apparent in the 
transform array F(n). In other words, transform array F(n) restricted to positive 
20 frequencies may be a poor approximation to the positive frequency image. Thus, the 
frequency location at which the magnitude fimction |F(n)| attains a maximum, when 
considered over positive frequencies, may only be a crude initial approximation to the 
tone frequency fj. 

Figure 12 is a plot of the magnitude of transform F(n) corresponding to a typical 
25 windowed input signal x(n), similar to that described above with reference to Figure 3. 
Note that the transform F(n) has a symmetry given by F(n) = F(k+N) for any integer k. In 
particular, F(-k)=F(N-k). Thus, frequency bin numbers between N/2 and N may be 
interpreted as negative frequencies. 
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In step 206, method may identify a frequency location proximate to an amplitude 
peak in the frequency transforai, wherein the amplitude peak may correspond to the tone. 
For example, in step 206 the method may determine the maximum amplitude peak in the 
input signal, which is presumed to be the peak of the tone. In one embodiment, the 
5 method may scan the DFT magnitude values |F(n)| over the range of positive frequency 
bins to determine the potential local positive maximxmi amplitude, e.g., to determine the 
bin index k which achieves the maximum magnitude. Where multiple tones are present 
in the input signal, the method may find multiple amplitude peaks and select one for 
subsequent processing. 

10 In step 208 the method may select two or more frequency bins proximate to the 

identified frequency location in the frequency transform. Thus the method may select 
two or more bins proximate to (or around) the positive maximum determined in 206. The 
method preferably selects frequency bins located on either side of the frequency location 
of the amplitude peak. In other words, at least one of the two or more frequency bins is 

15 on each side of the frequency location, i.e., one or more bins on a first side of the 
frequency location of the amplitude peak and one or more bins on the other side of the 
frequency location of the amplitude peak. In one embodiment, in step 208 the method 
may select several bins k that are in the neighborhood of k^, wherein k^ is the integer 
bin index value k in the range from 0 to N/2 which maximizes the magnitude of F(n). 

20 As noted above, the maximum of magnitude fimction |F(n)| considered as a 

fimction of continuous frequency typically does not occur at the integer value k^, but 
should occur in the interval between k^ and the next highest amplitude sample 
frequency, kj Figure 13 illustrates a blowup of the positive frequency magnitude peak 
301 proximate to bin index value k^^. In this example, the method may select k^^-l and 

25 kn^+1 as the bins proximate to the frequency producing the positive maximum 
amplitude. In this embodiment, the method selects three frequency bins for processing. 
However, the method may select 2, 3, 4, or a greater number of bins, as desired. 
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In step 210 the method may detemiine a tone frequency value that minimizes a 
difference between two or more expressions, e.g., at least a first expression and a second 
expression. Stated another way, the method may involve selecting a tone frequency value 
that makes a plurality of expressions most nearly equal. 
5 Each of the plurahty of expressions may comprise a simi of one or more 

numerator terms divided by a sum of one or more denominator terms, wherein each of the 
plurahty of expressions includes a tone frequency variable, wherein each numerator term 
and each denominator term corresponds to one of the frequency bins, and wherein a ratio 
of each mmierator term and its corresponding denominator term represent a complex 

10 ampUtude of the tone at a respective bin. 

Thus, for example, the first expression may comprise a sum of two or more 
numerator terms divided by a sum of two or more denominator terms, wherein the first 
expression includes a tone frequency variable. In a similar manner, the second expression 
may comprise a sum of one or more of the numerator terms divided by a sum of one or 

15 more of the denominator terms, wherein the second expression includes the tone 
frequency variable. The nimiber of numerator terms and corresponding denominator 
terms in the first expression are preferably different than the number of numerator terms 
and corresponding denominator terms in the second expression. 

In each expression, each numerator term and its corresponding denominator term 

20 correspond to one of the frequency bins. In one embodiment, a ratio of each numerator 
term and its corresponding denominator term represent a complex amplitude of the tone 
at a respective bin. In each expression, the tone frequency variable may represent a 
correct tone frequency value of the tone. The first expression may be approximately 
equivalent to the second expression when the correct tone frequency value is used for the 

25 tone frequency variable in the first and second expressions. 

As one example of the expressions used, where two frequency bins are selected, 
the first expression may comprise a sum of two mmierator terms divided by a sum of the 
corresponding two denominator terms, wherein the two numerator terms and the two 
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denominator terms correspond to the two frequency bins. The second expression may 
comprise one nimierator term divided by one corresponding denominator term 
corresponding to one of the two frequency bins. 

As another example, where three frequency bins are selected, the first expression 
5 may comprise a sum of three numerator terms divided by a sum of the corresponding 
three denominator terms, wherein the three numerator terms and the three denominator 
terms correspond to the three frequency bins. The second expression may comprise a 
sum of two of the three numerator terms divided by a sum of two of the three 
corresponding denominator terms, corresponding to two of the three frequency bins. 

10 The step of determining a tone frequency value may comprise computing a 

plurality of differences between the first expression and the second expression for 
different respective tone frequency values of the tone frequency variable, and then 
selecting the tone frequency value that produces a smallest difference. The method may 
compute the plurality of differences by performing a Newton-Rhapson root finding 

15 method, as is known in the art. 

The expressions may be real expressions or complex expressions. Where the 
expressions are complex expressions, the method may involve minimizing a difference 
between an amplitude of the first complex expression and an ampUtude of the second 
complex expression. 

20 

In one embodiment, the first expression and the second expression have the form: 

k+M , V k^M-l V 

tXF{n)W'iJ„ -f^)-F\n)W{f„ +/)) (/„ - f^)-F' {n)W{f„ +/)) 

_ k+\ 

Jt+A// \ k-¥M-\/ \ 

it+i Jt+i 



25 wherein: 
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F(n) is the nth value of the single sided scaled FFT spectrum; and 
W represents a window function, wherein the window function is shifted by a 
value of the tone frequency variable fj. 



Introducing the function G(f): 



G(f) = 



k+M . V 

Y,(nn)W\f, -f)-F\nWWn + A ^ 



k+l 

k+M 



k+M- 



Z W(j„-ff-\w(j„+f)\\ 



+M-} V k+M/ > 

2{nn)W'(J„ -f)-F\n)W(f„+fX x -ff -\nf„ +f)\\ 

k+\ k+\ 



k+\ 

k+M i 



10 



This results in: 



G(/) = 0 



15 



This equation is a function of f, and this equation can be solved by applying 
Newton Raphson*s root finding algorithm to determine the tone frequency/. The theory 
behind the use of these expressions is discussed further below. 

When the frequency transform of the samples computed in 204 comprises 
generating a power spectrum of the samples, the first expression and the second 
expression have the form: 



20 



k+M 

E(li^(«)|xK(/„-/,)|) 



k+\ 



k+M 



ZKa -fi)\ 



k+\ 



k+M- 



k+\ 



k+M-\ 

z 

it+1 



•riYi —1 

I \w{f„-Lt 
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wherein: 

F(n) is the nth value of the single sided scaled FFT spectrum; and 
W represents a window function, wherein the window function is shifted by a 
value of the tone frequency variable fj. 

5 

Let: 

This results in: 

G(/) = 0 

10 

Again, Newton Raphson root finding may be performed to solve for the frequency 
of the tone. The theory behind the use of these expressions is discussed further below. 

The tone frequency value determined in step 210 may comprise the correct 
15 frequency of the tone. As used herein, the term "correct" includes correct or substantially 
correct. Li other words, the determined tone frequency value represents a good 
approximation of the correct tone frequency value. The method described herein operates 
to find a more exact approximation of the frequency of the tone, even in the presence of 
noise or interference from other tones (i.e., the positive and/or negative images from other 
20 tones), or the negative image of the respective tone being found. The method is also 
more computationally efficient than prior art techniques. 

Where the input signal comprises a plurality of tones, the method may operate to 
correctly find signal parameters of a first tone of the plurality of tones. Where it is 
25 desired to locate two or more of, or all of, the tones present in the input signal, steps 206 
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-210 may be iteratively performed for each of the pluraUty of tones to determine at least 
one signal parameter for each of the plurality of tones in the input signal. 

After the tone frequency value has been determined in step 210, in step 212 the 
method may then optionally compute one or more of the amplitude and phase of the tone 
using the determined tone frequency value. For example, wherein the FFT spectrum is 
available, once the / value is known, the amplitude and phase can be computed as: 

(/„ -f^)-F\n)W(f„ +/;.)) 

Where the FFT power spectrum is available, once the frequency is known, the 
amplitude can be computed as: 

SH«)|x|fF(/„-/.)|) 

A = _^±1 

In step 214 the method may then comprise storing the determined tone frequency 
value, as well as the amplitude and/or phase if these are computed, in a memory, or 
outputting the determined tone frequency value, amplitude and/or phase, e.g., on a 
display. 

Theory of Operation of the Expressions 

The following describes the derivation and theory behind the expressions used in 
the preferred embodiment. 
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As noted above, given a single tone signal in the time domain, the input signal 
comprising the tone can be expressed in following fomi: 

x(n) = 4. cos(2;^7i + q>.) (1) 
Where Ai is the amplitude, / is the normalized frequency, n is from 0 to N-l and 
5 (pi is the initial phase. Here the subscript "i" is used to represent a single tone in what 
could be a multitone signal. 

If the FFT is computed based on the windowed version of x(n) as performed in 
step 204, then in the frequency domain, as noted above, each FFT bin can be computed as 
10 follows: 

Fin) = A,e-''^W{f„ +f^) + A,e'^'W(f„ (2) 

In the above equation, the first term of the equation represents the negative image 
of the tone and the second term represents the positive image of the tone. 

Each equation can then be converted to a different form. For example, given the 
following: 

C = V' (3) 

Then 

F(n) = CX/„ + + CWif, - f,) (4) 

Here, F(n) is the «th value of the single sided scaled FFT spectrum. 
The windowed frequency response can be easily computed numerically once the 
window coefficients are known. 

25 Estimate Amplitude/Frequency/Phase based on FFT spectrum 

If the FFT spectrum is available, each particular single tone can be roughly 
determined by finding the local maximum in the sense of amplitude. Here assume 



15 



20 
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F(k+1), F(k+2),.. F(k+M) are the M bins around the single tone. This resuhs in a set of 
nonlinear equations of : 



Fin) = C'W(f„ +f^) + CWif„ -/.) (5) 
5 Where n = k+1, k+2, . . . k+M; 

From this point, the following equation can be derived: 

^ F(n)W'(f„-f,)-F\rtWiL+f>) 
\nL-f<)f-\W{f„+f>)f 

For each n = k+1 , . . .k+M. 

10 Furthermore, to reduce the numerical fluctuation, the equation (7) may be 

written as follows: 



Z{Fin)W\f„ -fd-F\n)W(J„ +/)) - fi)-F\n)W{fn +/)) 

J+1 _ k+\ ^ 

k+M( ^ \ k+M-U \ 

jfc+i jfc+i 

15 

This equation is a function ofy; , and this equation may be solved in step 210 by 
applying Newton Raphson's root finding algorithm. 

Define G(f): 
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G(/) = 



k+M- 



k+M-\ 



jfc+l 



Jfc+l 

+A/-L J *+A// > 



(8) 



Now by solving the equation: 

G(f) = 0 (9) 
using Newton Raphson's root finding algorithm (or other suitable root find 
methods), the tone fi-equency can be determined 

Once the fi value is known, the amplitude and phase can be computed as follows: 



k+M J 



2{nn)W\f„ -f^)-F\nWiL +/>)) 
1 0 A,e'^' =C = 



Tinfn-fi)f-\nL+fif) 

Jfc+l 



(10) 



Estimate Amplitude/Frequency based on FFT power spectrum 

If the FFT power spectrum is available, then the amplitude and fi-equency of the 
15 tone may be estimated based on the FFT power spectrum. Since the phase information is 
not available in the power spectrum, the effect of the negative fi-equency cannot be taken 
into account. The equation for computing the fi-equency are as follows: 

k+M 

Z^F(n)\x\Wif„-f,)\) 

ifc+1 

k+M 

Jfc+l 



Jfc+Af-L 

Zr(«)IH^(/»-//)l) 



k+i 



k+M-\ 

ifc+1 



t-m — 1 



(11) 
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Let equation (12) be: 



/t+A/ k+M-l k+M'l k+M 

G(f)=^(\Fin)\x\W(f„-f)\)x X \W(J„-ff- I(iF(«)|x|Pf'(/„-/)|)xX|f>^(/„-/)f 

k+\ k+\ *+l 

5 This results in equation (13): 
G(/) = 0 

Again, the Newton-Raphson root finding method may be applied to solve for the 
frequency of the tone. 

10 Once the frequency is known, the amplitude can be computed by equation (13) as: 



4 = 



k+M 

Unn)\x\W(f,-f,)\) 



k+l 



k+M 
k+l 



(14) 



In one embodiment, the power spectrum may be averaged before applying this 
15 tone estimation. The measurement is more stable when using the averaged power 
spectrum. 



Thus, the above single tone estimation method may be used to find and 
characterize a tone in a received signal. Additionally, the method may be used to detect 
20 and characterize a plurality of different tones present within a signal by applying the 
method iteratively to find each respective tone within the input signal. As described 
above, in one embodiment, as each tone is determined or found, the determined tone may 
be subtracted out from the input signal to produce a modified input signal, and the 
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method may be then performed on the modified input signal to locate the next tone, and 
so on. 

It is noted that although the single tone estimation method described above (or 
variants thereof) is preferably used in the present invention, other single tone estimation 
5 methods are also contemplated, such as that described in U.S. Patent Application Serial 
No. 09/753,164, titled "System and Method for Estimating Tones in an Input Signal", 
which is incorporated by reference above. 

Applications 

10 Embodiments of the present invention may be used in various applications. In 

general, embodiments of the present invention may be used in any system where it is 
desired to detect tones, e.g., sinusoidal tones, present in a signal, e.g., where it is desired 
to detect the precise frequency, amplitude and/or phase of the tones present in the signal. 
For example, an embodiment of the present invention may be used in a DTMF (Dual 

15 Tone Multi-Frequency) system for detecting tones present in a signal, such as a signal 
generated by a keypad of a telephone. Embodiments of the present invention are also 
contemplated for use in applications involving sonar, radar (e.g. Doppler radar), 
frequency-shift keying applications, mechanical systems analysis, etc. For example, the 
reflections generated by multiple moving objects in response to a radar pulse have 

20 distinct frequencies dependent on their radial velocities with respect to the radar station. 
Thus, the frequencies of the reflections are usable for tracking the multiple moving 
objects. In another example, a mechanical system excited with a physical stimulus (e.g. 
an impulse) may manifest vibrations at one or more frequencies. The frequency, 
amplitude and/or phase of these vibrations may provide information to a system analyst 

25 about the nature of flaws in the mechanical system. Embodiments of the present 
invention may be used in a wide variety of applications, i.e. in any application where it is 
desirable to identify one or more tones present in an input signal. The above-mentioned 
appHcations are merely representative examples. 



Atty.Dkt.No.: 5150-62101 



Page 42 



Meyertons, Hood, Kivlin, Kowert & Goetzel PC 



Example Graphical Program for Single Tone Detection 

Figures 14A -18B comprise screen shots of a portion of a graphical program 
which implements one embodiment of the invention. Figures 14A -18B comprise screen 
5 shots of a portion of a graphical program or Virtual Instrument (VI) written in the 
Lab VIEW graphical programming language. The graphical program is hierarchical in 
nature, where Figures 14A and 14B comprise the top level front panel and block diagram, 
and figures I5A/B, 16A/B, 17A/B, and 18A/B are various sub programs or sub- Vis of the 
program. Various other sub programs or sub- Vis of this graphical program that are not 
10 necessary to an understanding of this embodiment of the invention are not included for 
convenience. 

Figures 14A and 14B illustrate a top level front panel and block diagram, 
respectively, of the graphical program titled "estimate single tone". This graphical 
program includes a block titled "single tone info" ("single tone information"). 
15 Figures 15A and 15B illustrate a front panel and block diagram, respectively, of 

the graphical program block titled "single tone information" contained in the diagram of 
Figure 7B. This graphical program includes a block titled "tone est." ("tone estimation"). 

Figures 16A and 16B illustrate a front panel and block diagram, respectively, of 
the graphical program block titled "tone estimation" contained in the diagram of Figure 
20 SB. This graphical program includes two blocks labeled "G(f)" and "G'(f)", which are 
used to compute the signal parameters characterizing a tone, as described above. 

Figures 17A and 17B illustrate a front panel and block diagram, respectively, of 
the graphical program block titled "G(f)". 

Figures ISA and 18B illustrate a front panel and block diagram, respectively, of 
25 the graphical program block titled "G'(f)". 

Thus, various embodiments of the invention provide means for identifying and 
estimating signal parameters characterizing multiple tones in a received signal. The 
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systems and methods described herein provide for improved estimation of the signal 
parameters for each of the multiple tones in the signal by iteratively removing the effects 
of other tones from each tone and applying a single tone estimation technique to estimate 
the signal parameters of the tone. The method may iterate over the tones multiple times 
5 to increase the accuracy of the tone parameter estimates, producing a set of final signal 
parameter estimates for each tone. 

Although the system and method of the present invention is described in 
connection with several embodiments, it is not intended to be limited to the specific 
10 forms set forth herein, but on the contrary, it is intended to cover such altematives, 
modifications, and equivalents, as can be reasonably included within the spirit and scope 
of the invention as defined by the appended claims. 
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